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Abstract — In this paper, we propose an optimal relay power 
allocation of an AmpIify-and-Forward relay networks with non- 
linear power amplifiers. Based on Bussgang Linearization The- 
ory, we depict the non-linear amplifying process into a linear 
system, which lets analyzing system performance easier. To obtain 
spatial diversity, we design a complete practical framework of a 
non-linear distortion aware receiver. Consider a total relay power 
constraint, we propose an optimal power allocation scheme to 
maximum the receiver signal-to-noise ratio. Simulation results 
show that proposed optimal relay power allocation indeed can 
improve the system capacity and resist the non-linear distortion. 
It is also verified that the proposed transmission scheme outper- 
forms other transmission schemes without considering non-linear 
distortion. 

I. Introduction 

Relay-assisted communication is a promising strategy that 
exploits the spatial diversity available among a collection of 
distributed single antenna terminals for both centralized and 
decentralized wireless networks. In most relaying networks, 
a two-stage relaying strategy is used lfn il22l . In the first 
stage, a source transmits and all relays listen; in the second 
stage, relays cooperatively forward the source symbols to 
the destination. Amplify-and-Forward (AF) and Decode-and- 
Forward (DF) are two most common relaying schemes HlJ. 
DF scheme always using cyclic redundancy check (CRC) will 
cause interruptions when the relay detects errors from the 
received message. If the length of source data block is signifi- 
cantly large, DF will incur unaffordable decoding latency f2l|. 
AF scheme only amplify the received signal without any other 
processes so that little latency and complexity are introduced. 
On the other hand, although AF incurs noise propagation 
which lowers the system performance, at high power regime 
AF based relay network also can achieve full diversity as DF 
based one |3 1. Therefore, AF based relay networks recently has 
attracted many attentions |4|-|6|. We also focus our attention 
on AF based relay networks. 

On the other hand. Orthogonal Frequency Division Multi- 
plexing (OFDM) has been absorbed into various wireless com- 
munication standards, such as Long Time Evolution (LTE), 
IEEE 802.16 (WiMAX) and 802.11. This is because OFDM 



has many benefits fit for nowadays high-speed data require- 
ments, e.g., high spectral efficiency, resiliency to multi-path 
distortion and Radio-Frequency interference [7|. As may be 
expected, OFDM technology has been applied into wireless 
relay networks to achieve more advantages |l8]-||9]. Cross- 
Layer also has been investigated in flSl-fSS). However, it is 
well-known that OFDM suffers from the high Peak-to- Average 
Power Ratio (PAPR) problem, which always makes the power 
amplifier (PA) work in its non-linear region and incurs non- 
linear distortion for input signal [7J. Consequently, PAPR 
also becomes a challenge to OFDM based relay networks, 
especially in AF based relay networks. The effect of non- 
linear PA for AF based relay network was investigated in [ 10 1 . 
And in |11|, two receiver techniques for nonlinear amplifier 
distortion compensation was proposed in an OFDM relay- 
assisted cooperative communication system. So far, there are 
little works involved in AF based relay networks with non- 
hnear PA. 

In this paper, we intend to investigate the impact of non- 
linear PA to the Nonlinear optimization of AF based relay 
networks. Particularly, the optimal relay power allocation is 
provided in AF based relay networks with non-linear PA. First, 
we introduce the Bussgang Linearization Theory which can 
model a non-linear process into an equivalent linear model. 
By this theory, we can obtain the average signal-to-noise 
ratio (SNR) of single source-relay-destination channel. At the 
destination, we also provide a practical non-linear distortion 
aware receiver, where Maximal-Ratio Combining (MRC) is 
employed to obtain spatial diversity. Besides, Optimal equiva- 
lent channel estimation through non-linear PA is also designed. 
Consequently, the equivalent received SNR is derived. To 
maximize average received SNR, we set up an optimal relay 
power allocation problem with total relay power constraint. 
After carefully analyzing, we found that if all relay power fall 
into a valid interval, the optimization problem is a concave 
function of relay power Therefore, we finally obtain the 
optimal relay power allocation strategy based on Lagrange 
Multiplier Method. Through simulations, we can see that the 
proposed power allocation outperforms other power allocation 



Fig. 1. Time slots of signal transmission 



Fig. 2. ISLA model 



Strategies where non-linear distortion caused by non-linear 
power amplifier is not considered. 

II. System Model 

In this section, we introduce all involved models in our 
following analysis. First, we need identify the signal model of 
AF relay networks. 

A. Signal Model of AF Relay Network 

We consider a wireless network with N relays, one source 
and one destination. Every node has a single antenna that can 
not transmit and receive simultaneously. Denote the channel 
coefficient from the source to the ith relay as fi and the 
channel coefficient from the ith relay to the destination as gi. 
Assume that all /; and gi are independent complex Gaussian 
random variables with zero-mean and variance 51^ and Sf^, 
respectively [4J. Note that we suppose there is no direct 
channel between the source node and the destination node. We 
further assume a block fading channel model, where channel 
gain stays constant during a time block and changes from 
block to block. It is assumed that the instantaneous channel 
is unknown to the transmitting node but perfectly known 
at receiving node [l]. Moreover, all relays are synchronized 
during relaying phase. The impact of synchronization error 
between relays is also beyond the scope of our discussion [3]. 

The signal transmission is divided into two phases (See 
FiglHi. During the first phase, the source broadcasts its signal 
s to all relays, where s could be a OFDM symbol ifTOll . Then 
the zth relay receives 



(1) 



where Ps is the source transmit power and Ui is the receiver 
noise which follows the complex Gaussian distribution with 
^ C7V(0,iVo)- Assume that there is E{ss*}=l. After that, 
the relay must amplify the received signal as 



PsSl+No 



(2) 



So, there is — Pi- In previous works with a linear 

power amplifier, the ith relay would forward Xi to the destina- 
tion during the ith time-slot of the second phase (See Fig. 1). 
That is to say the linear amplifier seems to be transparent for 



the transmitted signal. However, in our consideration, could 
access the non-linear region of the front-end power amplifier 
with a non-vanishing probability. Model the behavior of power 
amplifier as a non-linear function Fa{x), therefore, the signal 
received by the destination during the ith time-slot of the 
second phase is given by 



Vi = giFA^Xi) + • 



(3) 



where Vi is the Gaussian noise with vt ^ CJ\f{0,No). Note 
that we merge the non-linear effect of the source amplifier into 
relay amplifiers [lOJ, so that the source amplifier is assumed to 
be linear. When the destination collects all these N forwarded 
signals, the transmitted signal from the source node could be 
detected. The receiver algorithm is discussed in section III. 

B. Non-linear Amplifier Model 

We assume the power amplifier (PA) of a relay is mem- 
oryless. Thus the frequency response of a PA is flat for all 
feed frequency. A memoryless amplifier can be determined 
by Amplitude to Amplitude (AM/AM) and Amplitude to 
Phase (AM/PM) characteristics |13|. Denote AM/ AM trans- 
form distortion as and AM/PM transform distortion 
as Ap(|a;|), then the output signal of PA is 



FAix) = A„,i\x\)e^^''''3(x)+AA\x\)) 



(4) 



where arg(a:) means the angle of complex signal x. Various 
amplifiers are depicted in llT3l . In this paper, we choose the 
ideal soft-limiter amplifier (ISLA) model as an example to 
reveal our methodology. More practical and comphcated PA 
models are considered in our ongoing work. 
ISLA Model: QS 



AM) - { 

Ap{\x\)^0 



c, 



\x\ > Asat 



(5) 



where Asat is the input saturation amplitude, Aq is the output 
saturation amplitude and there is Asat — {Aq — c)/b. Fig. |2] 
shows the ISLA model. Without lose of generality, we set 

c = 0, 6 = 1 and Asat = Aq. 



C. Bussgang Linearization Theory 

Bussgang Linearization Theory lfT4l says that the output of 
a non-Hnear power amplifier can be expressed with a scaled 
version of input signal and an additive uncorrected noise term 
if the input signal is a Gaussian variable: 



Fa (x) — ax + d 



(6) 



where a is the linear scale factor (LSF), d is the non- 
Hnear distortion (NLD) and there has to be E{x*d} = 0. 
If the number of sub-carriers is large enough, d becomes 
Gaussian and d ^ CJ\f{0,5'^) ifTSll . Thus, take into account 
this condition, it is easy to derive following two corollaries. 
Corollary 1: 

E{x*Fa{x)} 



E{\x\^} 



Corollary 2: 



6l = E{\Fa{x)\^} - aE{xFA{xy} 



(7) 



(8) 



Observe the power amplifier model, we just need consider 
AM/AM transform in following analysis. Moreover, by (|2) we 
can see that the Xi follows the Gaussian distribution so that 
Bussgang Linearization Theory is available in our model. 

D. Single Channel SNR Model 

Apply the linear model of (|6]l, then the receiver signal of 
the «th time-slot at the destination can be expressed as 



Vi = 



Nn 



(9) 



where hi — aifigi is the equivalent channel coefficient with 
h ^ CAA(0, |ai 1252.52^) and = ai^J'^^^^^g.n^ + 
di + Vi is the equivalent receiver noise, i.e., Wi ^ 
CM (0, '"pf/'l^^" + + iVo) ■ As a result, we can define 
the average signal-noise ration (SNR) of the ith received signal 
as 



(10) 



where pi 



V - 



iQfiP — 1 and 5'fi ■ — 0, then p 



If the power amplifier is linear, i.e., 
No/Pi and ( fTOb becomes 
coincident with the SNR expression in 15] ■ According to 
Corollary 1 and 2, there is 

_ P,iE{\FA{x,)\^} + No) 



\E{x*Fa{x,)W 



- 1 



(11) 



In the Appendix, we obtain the exact expression of pi based 
on ISLA model. Therefore, we have: 



1 - 



,V7lil 

where Erfc(i) — ^ 
ized relay power and p 



^Erfc 



(12) 
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_ No 



^dx, ^i = ^^5— is the normal- 
is the normalized noise power 



III. Optimal Relay Power Allocation 

After all TV time-slots of phase 2, the receiver at the 
destination intend to exploit {yi, j/at} to detect the signal s 
from the source. Obviously, if we need obtain high power- 
efficiency, the effect of NLD should be considered in the 
design of receiver 

A. Receiver Design 

A NLD-aware receiver was designed in ifTTI , where 
Maximal-Ratio Combining (MRC) is employed to obtain spa- 
tial diversity. However, how to compute ai and ^ and how to 
perform MRC with received signals are not provided therein. 
In this paper, we also adopt the NLD-MRC idea. Different 
from |11|, we provide a complete procedure of NLD-MRC 
receiver and a practical solution to estimate the parameters 
of the linear model (a and 5'jj). Like Il3|-|l5l, we assume 
the statistical channel state information and power allocation 
strategy are known by the receiver all the time. Further more, 
we also assume the receiver has the knowledge of the non- 
linear power amplifiers equipped at relays. 

Observe (|9]l, to perform NLD-MRC receiving, the receiver 
needs know each equivalent channel coefficient hi. We intend 
to use the pilot symbol Sp from the source node to estimate 
{hi}. Denote the received pilot signal as yi{sp) during the 
ith time-slot. Because the equivalent noise Wi is Gaussian, the 
optimal estimation of hi is [16] 



(13) 



where Sl^ = "^^^^ + J^, + No and Si = la^lH^Sf,. 
Obviously, we need firstly estimate {a^} and {6^ ^} to obtain 
the equivalent channels {hi} at the NLD-MRC receiver For 
this purpose, we provide a following receiving algorithm: 



1) Step 1: Use the results of Appendix to compute 
E{x*FAix^)} and E{\FAix^)\^}. 

2) Step 2: According to two corollaries in Section II, obtain 
the linear model parameters {a,} and {(5^ j}. 

3) Step 3: By ( fT3b . estimate the equivalent channel coeffi- 
cients {hi} 

4) Step 4: Detect the transmitted signal with rule: 



arg mm ■ 



N 



1=1 Z^i=l 



In following context, we consider an ideal case where hi = 
hi to analyze the optimal system performance. In term of the 
NLD-MRC receiver, the average received SNR is 



N 

E 

j=i 



r,; 



At this point, both selection combining and equal-gain comb- 
ing could also be employed herein. To achieve best perfor- 
mance, we only consider MRC in this paper. 




B. Relay Power Allocation Problem 

Assume all power amplifiers at relays have the same model 
and parameters. Observe (ITOt . we can see that Pg and {Pi} 
are the variable factors to optimize SNR. As stated in flS], 
it is difficult to optimally allocate Pg and {Pi} simulta- 
neously. We also consider a relaxed problem to one with 
a fixed pre-determined P,. For a relay, the actual average 
transmit power is E{\FA{xi)\'^}, which is a function of Pi. 
Assume the maximum total relay power as P,., we have 
E^i ^Lt7» (l - e"^) < Pr by Define the normal- 
ized maximum total relay power constraint 7^. = Pr/A^^f.. 
Therefore, the optimal power allocation problem is 



max {Ft } 

N 

8.1.^7,; (1 - er~ ) < Jr 



i=l 



(14) 



0<7„ i^l,...,N 



Observe Fig. [3] although F; is not a concave function of 
7i, Ti could achieves the unique maximum value if = 0. 
Then we can derive the optimal 7^ without power constraint 
is the solution of 



1 



2fi 



-Erfc 



1 



(15) 



By ( fTSl ). we can also confirm that there exists an unique 
maximum value. We denote the solution of ( flSl l as 7op((/i) 
which is independent of channel state information. Through 
Fig. (|3]l, we can see that if 7^ becomes larger than jopt, Tj 
decreases. To save power, we have no need to allocate power 
larger than 'jopt to lower the receiver SNR. Therefore, we just 
need consider < 7,; < 7opt as an efficient interval lfT2l . 
In addition, we plot '^opt versus /i curve in Fig. |4] It shows 
that 7opt is a monotonic function of /i. Therefore, we could 
set up a lookup table to show the value of '^opt for different 
normalized noise power Furthermore, we found that 
if 7i falls into the interval < 7^ < 'jopt- 



< 



Off 

Thus, the Hessian 



matrix 



V^Ft = diag 



N 



dll ' 97I ^7^ 



(16) 



is a negative semidefinite matrix in the efficient interval. It is 
obvious that Ff is a concave function [17] of (71, 72---, 7jv) 
if < 7i < 7opt, i = I, ...jN. Therefore, we can rewrite the 
optimal power allocation problem as 



N 



S.t. 



7i (1 - e"^) < 7r 

< 7i < 7opt, i = l,- 



(17) 



First, relax the optimization problem with only total relay 
power constraint. Construct the Lagrange function 



N / N 

L = Y.^^ (7.) + A ^7.(l-e-* 



(18) 



where A > is the Lagrange factor Let ^ = 0, we obtain 
(ISTT i. To simplify the expression, we denote the equation of 
(ISTT i as A = ^{ji,6gj^,dfj). Hence we can express 



(19) 



where $ ^(0 is the inverse function of $(.), A is a constant 
to meet the total relay power constraint and 





X 



if a; < 
else 



Recall Fig. [3] F, is a monotonically increasing and concave 
function of 7i if < 7^ < 7opt. Therefore, the optimal solution 
of ([TtT i must be on the boundary |5|[17|. As a result, the 
optimal relay power allocation of (flTt is 



lopt 



(20) 



[Sf^No + Pi{Ps6l^ + No)f 2'yf 



V^JJ 



^Erfc 

7! 



(21) 




SNR (dB) 

Fig. 5. Capacity of the symmetric relay network in four cases 

where 

f if a; < 

(^)o ~ { ^ < X < y 

[ y if x> y 

It is obvious that the optimal power allocation solution can 
be considered as an extended water- filling process |l5], where 
there is a lid and bottom for each vessel. The lid is generated 

by Jopt{^J.)■ 

The destination can estimate and collect all statistical chan- 
nel state information of the relay network. Since the model 
of power amplifier at relays is also known by the destination, 
the destination can obtain the optimal relay power allocation 
by (I20I 1. Through feedback link, the optimal allocation can be 
informed to all relays. As the optimal power allocation is based 
on the statistical channel state information, computing for 
optimal power and feedback link are not always be activated. 
If one relay is assigned with zero power, this relay will not 
forward signal during the phase 2. To enhance the spectrum 
efficiency, the destination could adaptively program the time- 
plots of phase 2 for relays with non-zero assigned power 

IV. Simulation Results 

In this section, we show simulation results to verify our 
analysis and to show performances of optimal power alloca- 
tion. We simulate an AF network with a source, a destination 
and four relay nodes {N = 4). Without lose of generality, 
assign source power as P5 = 1 and noise variance as 
N„ = -15 dB Cni. Set A^at = h then fj. -15 dB and 
Through (flSl l. we have 'jopt = —2.8507 dB. Define the relay 
network SNR as Pr/N^. Four transmission schemes with the 
same power constraint are compared: 



8 




1 1 1 1 1 1 1 

5 10 15 20 25 30 

SNR (dB) 

Fig. 6. Capacity of the non-symmetric relay network in four cases 



• With linear PA, traditional water-filling ||5] is used (LPA- 
TWF). 

• With Non-linear PA, traditional water-filling |p] is used 
(NLPA-TWF). 

• With Non-linear PA, equal power allocation {Pi = Pr/N) 
is used (NLPA-EPA). 

• With Non-linear PA, proposed power allocation is used 
(NLPA-Proposed). 

In Fig. |5] we simulate a symmetric relay network, where 
each channel variance is unit, i.e., ^ CAf{0., 1) and Sf^ ^ 
CAf{{), 1). Obviously, in this case, NLPA-EPA and NLPA- 
TWF converge into an identical case, which is also verified 
by Fig. |5] As there is no NLD in LPA-TWF transmission 
scheme, it achieves the best performance in all four transmis- 
sion schemes. So, LPA-TWF is an ideal transmission scheme 
for other three schemes. To overcome NLD, a transmission 
scheme with non-linear PA should approach the performance 
of LPA-TWF as close as possible. The more close, the more 
ability of resisting NLD. From 5 dB to 15 dB, all transmission 
schemes achieve the same capacity. The reason is that the relay 
power is operated in the linear region of PA and the NLD 
derived by Bussgang theory is very slight, i.e., pi « No/Pi. 
If 15 dB < SNR < 18dB, all other three schemes except 
LPA-TWF have a slower and slower capacity growing rate as 
SNR increases, which is because NLD starts to violate the 
transmitted signal from relays and the relay power gain is 
still larger than the harm of NLD. The higher relay transmit 
power is, the heavier NLD is caused. As the SNR further 
increases, if SNR > ISdB, the capacity of LPA-TWF still 
increases and NLPA-Proposed power allocation keeps a fixed 



capacity, but the capacities of NLPA-TWF and NLPA-EPA 
become to decrease. The reason is that each normalized relay 
power has exceeded the ^opt = —2.8507 dB and relay power 
gain becomes smaller than the harm of NLD. Therefore, it 
is verified that our proposed indeed overcomes part of the 
NLD and outperforms other power allocation schemes without 
considering NLD. 

We also consider a non-symmetric relay network to verify 
our proposed power allocation in Fig. |6] Set each channel 
variance is randomly selected from the interval (0.5, 1.5). We 
average our simulation results over 1000 channel variance 
realizations. Compare Fig. |5] and Fig. |6] we found that the 
only difference is that in both low SNR regime and high SNR 
regime NLPA-TWF outperforms NLPA-EPA. In other words, 
the advantage of water-filling algorithm based on statistical 
channel state information lies in the differences between 
channel variances. Through entire SNR region, the proposed 
power allocation achieves the best transmission capacity in 
contrast to NLPA-TWF and NLPA-EPA and there is also a 
capacity gap between NLPA-proposed and NLPA-EPA in low 
power regime. That is to say the proposed power allocation 
has a better performance in non-symmetric relay networks. 
In summary, the proposed optimal relay power allocation can 
improve the system capacity and resist the NLD incurred by 
non-Unear PA. 

V. Conclusion 

An optimal relay power allocation is proposed in this paper 
to maximize the receiver SNR and capacity. With the bussgang 
Unearization theory, we convert the non-linear signal incurred 
by non-linear power amplifier into an equivalent linear model. 
A complete framework of non-linear distortion aware receiver 
which can perform maximal-ratio combining is proposed to 
achieve spatial diversity. After that, we proposed a optimal 
power allocation scheme with total relay power constraint. 
The proposed power allocation outperforms all other power 
allocation schemes in the presence of non-linear power am- 
plifier and has the ability of resisting the non-linear distortion 
caused by non-linear power amplifier. Through simulations, we 
also found that proposed power allocation can achieve better 
performance in non-symmetric relay networks. 

Appendix 

By (|2]), Xi is also Gaussian, i.e., Xi ^ CJ\f{0, Pi). Then there 
is i^llxip} = Pi. Obviously, the amplitude of Xi follows the 
Rayleigh distribution and with probability density function 

fi\x^\) = -^^^P ( — 

. According to (|4]i, there are 

E{x*FAix.O} - 

/■+°° (22) 
/ |x,|A„(|a;,|)e^"^-(l-'l)/(|x,|)d|x,| 
Jo 

and 

r+oo 

E{\FAix,)f}= / A^{\x,\ff{\x,\)d\x,\ (23) 
Jo 



Consider Ap(|xi|) = 0, then we can deduce 

E{x*Fa{x,)} ^ E{x,FAix,)*} 
Apply the model of ISLA, we have 

E{x*Fa{x,)} = 



P 1 



exp 




fErfc 



A.. 



/P 



^^{|^a(a:.)l'} = ^'^ 1-cxp - 



P. 



(24) 



(25) 



(26) 
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